NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


Rear  Admiral  J.  j.  Ekelund  Jack  R.  Borsting 

Superintendent  Provost 


The  work  was  prepared  by: 


^D/P.  GaverV^Professor 
Department  of  Operations  Research 


P.  A.  Jacc±>£,  Associate  Professor 
Department  of  Operations  Research 


I 


Reviewed  by: 


Michael  G.  Sovereign, 
Department  of  Operations  Resei 


Released  by: 


William  M.  Tolies 
Dean  cf  Research 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Fnteret f) 


t 


1 


IX  aUTHf»H/.L 

%Y  p-/c 


— rt 


I  NPS55-  W7 
j  TiTiiC  [■■urTTitrrrrfrr- — 


REPORT  DOCUMENTATION  PAGE 


r 


2.  GOVT  ACCESSION  NO 

ktMEM. 


TITI  S  (■llf  71111111111 

On  Combinations  of  Random  Loads 


COVERED 


Gaver 


P.  A.^Jacobs^ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


6.  CONTRACT  OR  GRANT  NUMBER^; 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School  j 
Monterey,  CA  93940 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
ArtEA  A  WORK  UNIT  NUMBERS 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA  93940 


14.  MONITORING  AGENCY  NAME  ft  ADDRESS (ft  different  from  Controlling  Office) 


16.  DISTRIBUTION  STATEMENT  (of  thle 


15.  SECURITY  CLASS,  (of  thlo  rawer*; 

Unclassified 


ISa.  oeclassification/oowngrading 
SCHEDULE 


Approved  for  Public  Release;  Distribution  Unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ahatracl  entered  In  Block  20.  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  rereree  aide  If  neceeeery  end  Identify  by  block  ntmber) 

Random  Load  Combinations, 

First  Passage  Times, 

Asymptotic  Properties, 

Maxima  of  Correlated  Random  Variables 


p.  ABSTRACT  (Continue  an  rereree  aide  II  neceeeery  end  Identify  by  block  number) 

S-  cuctures  are  subject  to  changing  loads  from  various  sources.  In  many  in¬ 
stances  these  loads  fluctuate  in  time  in  an  apparently  random  fashion.  Certair| 
loads  vary  rather  slowly  (called  constant  loads) ;  other  loads  occur  more 
nearly  as  impulses  (shock  loads) .  Suppose  that  the  stress  put  on  the  struc¬ 
ture  by  various  loads  acting  simultaneously  can  be  expressed  as  a  linear 
combination  of  the  load  magnitudes.  In  this  paper  certain  simple  but  somewhal| 
realistic  probabilistic  load  models  are  given  and  the  resulting  probabilistic 


DD  , 


FORM 
JAM  73 


1473 


EDITION  OF  1  NOV  66  IS  OBSOLETE 
s/N  0102-014*  6601  . 


OVER 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fW*>«i  Bpretnfered) 


UNCLASSIFIED 


‘  WCty***  ^^iK^'-'vxv^r,'  - 


ON  COMBINATIONS  OF  RANDOM  LOADS 


by 
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and 
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Monterey,  CA  93940 


The  integrity  of  a  great  many  physical  structures  is 


potentially  threatened  by  combinations  of  physical  loads  of  vary¬ 
ing  magnitudes  from  various  sources.  We  are  thinking  of  structures 
such  as  buildings  (for  instance  those  that  house  or  contain 
nuclear  power  plant  elements) ;  aircraft  and  spacecraft;  electrical 
transmission  networks;  bridges,  piers,  and  dams;  and  offshore 
oil  drilling  rigs  that  experience  loads  from  wind,  snow  and  ice, 
tides,  earthquakes,  and  so  forth.  In  many  instances  the  total 
load  or  stress  experienced  by  a  structure  varies  in  time  in  an 
apparently  random  fashion.  Certain  load  components  vary  rather 
slowly;  for  example,  that  component  resulting  from  snow  and  ice 
accumulation;  others  occur  more  nearly  as  impulses,  such 
as  those  associated  with  winds  or  earthquakes.  The  problem  is 
to  design  structures  to  withstand  the  superposition  of  loads 
from  many  sources  with  at  least  an  approximately  understood 
(high)  probability.  In  engineering  terms  we  wish  to  work  towards 


C 

V4 


1 


developing  a  rational  safety  factor  technique  for  designing 
structures  to  withstand-  the  combination  of  loads  anticipated. 

The  purpose  of  this  paper  -is  -to  describe  and  investigate  certain 
simple  but  somewhat  realistic  probabilistic  load  models  for 
use  in  design,  and  perhaps  safety,  assessment  of  structures. 

In  this  paper  we  confine  attention  to  the  superposition 
of  just  two  load  types:  shock  loads ,  and  constant  loads.  For 
example,  wind  gusts,  flash  floods,  and  earthquakes  have  varying 
magnitudes  and  have  relatively  short  durations  in  comparison  to 
the  times  between  their  occurrences;  these  will  be  modeled  as 
instantaneous  shock  loads.  On  the  other  hand  snow,  ice,  or 
water  accumulation,  or  even  the  presence  of  slow-moving  vehicles 
or  furniture,  present  loads  that  remain  nearly  constant  in  time, 
occasionally  changing  to  new  levels;  these  will  be  modeled  as 
constant  loads  that  change  infrequently.  Throughout  this  invest! 
gation  it  will  be  assumed  that  the  effective  stress  exerted  by 
several  types  of  loads  acting  simultaneously  can  be  expressed 
as  a  linear  combination  (actually,  sum)  of  those  component  loads, 
the  loads  being  treated  as  stochastic  processes.  Somewhat 
similar  load  combination  models  have  been  studied  in  the  past; 
see  the  work  of  Bosshardt  [1975],  Larrabee  and  Cornell  [1978], 
McGuire  and  Cornell  [1974],  Pier  and  Cornell  [1973],  ,Weh  [1977], 
and  others.  Note  that  the  present  study  utilizes  the  notion  of 
stochastic  point  processes  as  a  modeling  tool;  see  Cox  and  Miller 
[1965]  for  an  introduction,  and  Lewis  (ed.)  [1972]  for  a  more 

extensive  treatment. 


Modeling  Assumptions  and  Problems 


Assume  that  the  successive  magnitudes  of  shock  loads  are 
independently  and  identically  distributed  (i.i.d.)  random  variables 
with  common  continuous  distribution  function  G(y);  G(y)  =  0 
for  y  £  0.  The  shock  loads*  appearance  is  regulated  by  a  Poisson 
process  with  rate  P.  Furthermore,  the  constant  loads  have  i.i.d. 
magnitudes  with  distribution  F(x) ;  F(x)  =  0  for  x  <  0.  Constant 
loads  change  at  moments  of  a  Poisson  process  with  rate  A . 

Finally,  the  shock  and  constant  load  processes  are,  for  the 
most  part,  taken  to  be  statistically  independent,  although 
this  assumption  can  be  relaxed  at  times  without  difficulty. 

Let  X(t)  (respectively  Y(t))  be  the  magnitude  of  the 
constant  (shock)  load  process  at  time  t;  (Y(t)  will  usually 
be  zero).  Put  Z(t)  =  X(t)  +  Y(t),  the  superposition  of  the  two 
loads  at  time  t,  and  M(t)  =  sups <  Z(s),  the  maximum  load 

combination  in  £0 , t] .  See  Figure  1.  Let  Tx=(inf  t^OrZft)  >  x} 

A  constant  load  process  I 

X(t)| 


SHOCK  LOAD  PROCESS 


LOAD  COMBINATION  PROCESS 


Fig.  1 

THE  LOAD  COMBINATION 


i.e.  the  first-passge  time  to  a  combined  stress  level  x.  Unless  fatigue 
or  some  other  cumulative  effect  is  in  operation  T  will  represent 
the  time  to  failure  of  a  structure  whose  strength  is  x,  and  that 
is  subjected  to  a  stress  history  (Z(t),  t  _>  0) . 

Section  2  is  concerned  with  the  description  of  the  distri¬ 
bution  of  the  maximum  process  and  the  distribution  of  T  .  The 

X 

first  step  is  to  obtain  Laplace  transforms  with  respect  to  time. 
Asymptotic  results  (x  -*■  «)  for  E[T  ]  and  the  distribution  of 

wv 

the  normalized  value  T*  =  T^fEfT^])  will  then  be  presented.  In 
particular  it  will  be  shown  that  the  distribution  of  T*  is 
approximately  a  unit  exponential  distribution  as  x  becomes  large 
and  so  does  E [T  ] . 

Section  3  is  devoted  to  a  study  of  the  asymptotic 
properties  of  the  "maximum  process,"  M(t),  as  t  -*■  These 
results  relate  to  those  of  Welsch  [1972]  and  O'Brien  [1974a, b]; 
they  are  seen  to  extend  the  classical  "extreme  value"  results 
of  Gumbel  [1958]  and  Gnedenko  [1943].  We  present  results  for 
the  asymptotic  distribution  of  M(t)  and  for  the  joint  distri¬ 
bution  of  the  first  and  second  maximum  load  combination  in 
[0,t]  as  t  ». 
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2.1.  Towards  the  Distribution  of  M(t) 


Consider  first  the  distribution  function  of  the  maximum 
process,  H^(t)  =  P{M(t)  _<  x} ,  for  x  >  0.  Its  Laplace  transform 
i?  available  immediately  by  noting  that  if  is  the  time  of 

the  first  change  in  the  magnitude  of  the  constant  load  process, 
then 


«>.,( t)  =  P{M(t)  <  x,  Tx  >  t}  +  P{M(t)  <_  x,  Tx  <  t} 


-At  *  - 

=  e-At  /  exp{-pt  G(x-y) }  F(dy) 

0 


t  -Av  x 

+  /  Ae  dv  /  exp{-yv  G(x-y)}  H  (t-v)  F(dy) 
0  0  x 


(2 


Next  take  Laplace  transforms  with  respect  to  t: 


=  l  e  Kt  Hx(t)dt  =  MvU)  +  AMv(?)  h„(5)  ; 


(2 


reversal  of  the  order  of  integration  provides  t:  vt 


M  (£)  “  /  U  +  A  +  yG(x-y)  ]  1  F(dy)  . 
x  0 


(2 
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(2.4) 


Hence, 

MxU) 

hx(?)  =  i  -  ajTTTT  * 

«?» 

I <  jeems  to  be  difficult  to  solve  (2.1)  (invert  (2.4))  in  any 
simple  convenient  form  for  any  interesting  choice  of  the  distri¬ 
butions  F  and  G.  However,  useful  information  can  still  be 
gleaned  from  (2.4).  First,  it  is  clear  that  if  T  is  the 
first-passage  time  to  combined  load  x,  so 


then. 


bet 


since 


F,  *  0 


Tx  =  inf(t  2  0:M(t)  >  x}  , 
P {'M ( t)  <  x}  =  P{Tx  >  t>, 

oo 


h(U  =  I  e~?t  p{t  >  tldt  > 
x  o  x 

in  (2.4:)  to  find  that 

M  (0) 

m(x)  eE[Tx]  =  hx(0)  =  r  -XM  W  * 


(2.5) 


(2.6) 


(2.7) 


We  now  record  some  expressions  for  the  mean  first-passage  time 
to  x  when  specific  distributions  for  shock  and  constant  load 
magnitudes  are  in  force. 
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Example  2.1.  Identical  exponentials:  P(x)  =  G(x) 
In  this  case 


x  > 


..  ,„s  1  /,  _-xx  yxe“x  ye”x  T  X  +  ye  x  "] 

Mx(°)  -  x  (1  -  e  )  -  —2 - —  in  |_-x  7  u  J 

and 

E[T  j  =  f  X  (l-e~x)  -  yxe~x  +•  -ye~x  Sin  [  (A+y),/(A+ye~x)  ] 
x  X2  j  1  +  &  x  +  ±  Hn[(X+y)/(X+  ye~x)  ] 

1  /ex\ 

~  y  \  T  J  as  x  -  -  . 


Example  2.2.  Different  exponentials: 


i 

i 

i 

k 

% 


F(x)  -  e-ax,  G(x)  =  e~kx  where  ^  =  2. 

In  this  case 

M  (0)  =  /X  [X  +  ye**k (x~y)  ]~X  ae~ay  dy 
x  0 

=  X-2  e”2kx  N (x) 

where 

N(»  =  X  [-1  +  e2kx]  +  2y[l-ekx] 

+  2y2  [-kx  +  £n[(X+y)(X  +  ye  kx)  ^]  ] . 


Thus 


E[T  ]  =  a”1  N(x)  [Ae2kx  -  N(x)]"1  ~  §- 
x 


(2.9) 


Example  2.3.  F(x)  =  e“x,  G(x)  =  e~x  .  In  this  case 


where 


Thus 


M  (0)  *  r  [1  '  e"XJ  +  £  e"x  N(x) 


N(x)  =  (e”x  -  1)  [(A  +  y)  (A  +  ue~x) ]-1 


E(Tx]  =  |ex  [1  -  e“x  +  ye'x  N(x))[l  -  yNfx)]''1  ~  I'v  e** 

(2.10) 


Example  2.4. 


F(x)  = 


G(x)  = 


x  <  aj/a  , 


x  >  a.1/*  , 


x  <  a. 


x  >  a. 


where  a  >  0  and  a.^,  a2  >  0.  In  this  case  as  x  » 


E(TX]  ~  [Aa1  +  ya2]"’1x0t. 


(2.11) 


Example  2.5. 


F(x)  = 


and 


G(x)  = 


c  a 

1 

.-3 


if 

if 

if 

if 


If  a  <  0 ,  then  as  x  -*■  » 


E[TX]  x“ 


If  0  <  a,  then  as  x  -*•  °° 


x  <  1  , 
x  1  , 

x  <  1  , 
x  >  1  . 


(2.12) 


The  next  result  is  a  limiting  result  for  the  first 
time  the  load  combination  process  exceeds  a  given  level  x. 

THEOREM  (2.1).  The  limiting  distribution  of  Tx  is  exponential 
in  the  sense  that 


lim 
x  -*■  x 


0 


PimU)-1^ 


>  t} 


where 


xrt  =  inf (t:F*G(t)  =  1} . 


Proof.  Put  p  =  pX-  and  rearrange  (2.3)  so  that  the 
denominator  of  (2.7)  is  in  the  form 

X  i 

1  -  XM  (0)  =  F(x)  +  /  pG(x-y)  [1  +  p§(x-y)  ]  F(dy)  rhjtx).  (2.13) 
x  0 

Now  apply  the  bounded  convergence  theorem  and  the  fact  that  P 
and  G  have  densities  to  show  that  m(x)  -*■  °°  as  x  -*■  Xg.  It 
next  follows  by  a  change  of  variables  in  (2.4)  that 


/  e~?t  P{m(x)~1  T  >  t}dt  =  m(x)-1  h  (y(x))  = 


.  -  M  (y(x)) 

/—V  \  —  A _ _ . 


(2.14) 


where  y(x)  =  £m(x)  .  By  the  bounded  convergence  theorem  applied 


to  (2.3), 


lim  M  (y(x))  =  y 

x**o  x  * 


(2.15) 


Again  from  (2.3)  we  have 


1  -  XMx(y(x)) 


-  F  (x)  +  /  [y(x)X  +  pG(x-y)  ]  [y  (x)  X  X  +  l+  pG(x-y) ]  F(dy) 
0 


Hence, 


m(x)  [1  -  AMx(y(x)  )  ] 

r  X  , 

=  X  /  tY(x)A“  +  1  +  pG(x-y)  ] -1  F (dy) 

x 

+  m(x)  £F(x)  +  /  pG(x-y)  [y(x) A~X  +  1  +  pG(x-y) ]“X  F(dy) ] 

Let 

h2(x)  =  F(x)  +  /  pG(x-y)  [y  (x)  A-1  +  1  +  pG(x-y)  ]~X  F(dy)  . 
then 

(x)-h2(x) 

x 

=/  F (dy)  PG(x-y) {- [1 +  pG(x-y) ]  +  fy (x) A_1 +  1 +  pG(x-y) ] } 

xUl  +  pG(x-y)  ]  [y  (x)  A_1  +  1  +  pG(x-y)]}"1. 
x 

‘l  FWy)  P°(x’y)  Y(x)A-1  ni+pG(x-y)]  [ylxiA-i  +  lt  G(x-y)]}-1. 
=  Y (x) A-1  k(x) 

where 

k<x)=/  pG(x-y){[l  +  pG(x-y)  ]{y(x)  A*1  +  I  +  pG(x-y)  J  }”yp(dy) . 
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By  the  bounded  convergence  theorem 


lim  k(x)  =0  . 
x  -*-xA 


Recall  that 


Hence, 


m(x)  = 


\  +  0(x) 
(x) 


hr(x)  -  h2(x)  Jh^OOT1  k(x) 


m(x) 


A*"1  +  0 (x)  hl(x) 


=  TT-olW  k(x) 


which  tends  to  zero  as  x  x^.  Therefore, 


lim  m(x)  [1  -  AM  (y(x) )  ]  =  £A  1  +  A-1  lim 
x  xQ  x  xQ 


1  + 


h2  (x)  -  h 


=  u+  D-r1. 


The  result  now  follows  from  (2.14),  (2.15)  and  the  unicity 
result. for  Laplace  transforms . 


3.  Asymptotic  Results  for  the  Maximum  Load  Combination 

In  this  section  we.  will  study  the  asymptotic  behavior  of 
the  maximum  load  combination  to  occur  during  time  interval  (0,t] 
as  t  ».  Most  results  concerning  maxima  of  non-Gaussian  random 
variables  are  for  variables  in  discrete  time.  Thus  we  will  first 
obtain  results  for  an  embedded  discrete  time  load  combination 
process . 


Let  Sn  be  the  time  of  the  nth  change  in  magnitude  of 
the  load  combination  process;  that  is. 


and  for  n  >  1 


=  inf {t  >  0:2 (t)  f  Z(t-)} 


=  inf {t  >  Sn-1:Z(t)  *  Z(t-)}, 


(3.1) 


The  change  in  the  load  combination  at  time  Sn  may  be  due  to 
either  a  change  in  the  constant  load  process  or  to  the  arrival 
of  a  shock  load.  Let  Zn  =  Z(Sn)  (respectively,  Xn  =  X(Sn) , 

and  =  Y(Sn))  be  the  magnitude  of  the  load  combination 

(respectively  constant  load  and  shock  load)  at  time  S  .  Put 

n 


M  =  max  (Z.)  . 
0  <  k  <  n  * 


(3.2) 


Note  that  the  times  (a^)  are  the  arrival  times  of  a  Poisson 
process  with  rate  A  +  u  and  are  independent  of  {zn>.  Further, 
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(3.3) 


P{Zq  £  x)  =  F(x) 


and  for  n;  >,  1 


where 


p{Zn  £  x)  =  pF(x)  +  qF*G(x) 


P  =  / '  *e 


r-Xt 


e  dt  =  t 


A  +  y 


and 


q  =  1-p.- 


3,1 .  Properties  of  An  Imbedded  Maximum  Process 


We  will  first  study  the  asymptotic  behavior  of  Mn  as 
n  -►  oo.  Note  that  the  random  variables  { Z^}  are  not  independent. 
In  fact  UX^Z^;  n  °  0,1,...}  is  a  discrete  time  Markov  pro¬ 
cess.  The  following  result  describes  the  dependence  of  the 
random  variables  (Zn>. 


PROPOSITION  (3.1) .  The  sequence  {Zn>  is  uniformly  mixing 
(cf.  Loynes  [1965,  p.  994]). 


Proof.  Let  A  and  B  be  events  such  that  A  -  f(Zg,...,Zn) 
and  B  =  g(Zn+In,  Zn+m+1,  ...).  Then 

|p  (A  ft  B)  -  P(A)  P  (B)  I 

=  |E[(P(B  |Xn)  -  P  (B)  )  ;  A]| 

=  qm|E[  /( ex  (dy)  -  F(dy))  P(B|Xn+m  =  y) ;  A] I 
n 

£  2qm  P (A) 

where  ex  denotes  the  Dirac  measure  concentrated  at  the  point  x. 
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Proof.  -  Note  that  for  any  c  >  0 

F(zi+l  >  >  <=>  > -J1  • 

The  .result  now  follows* 


In  order  to  obtain  the  limiting  distribution  of  Mn  we 
will  first  study  the  limiting  distribution  of  a  related  process 
which  is  defined  as  follows. 

Let 

1  if  Xn  -  0  - 
0  if  Yn  >  0  ; 

that  is,  Jn  =  1  if  the  nth  change  in  the  load  combination  process 
{Z (t) ;  t  2  Q}  is  due  to  a  change  in  the  constant  load  and  0  if 
it  is  due  to  the  arrival  of  a  shock  load. 


Let 


andr  s 


Tq  =  0 


xn  -  inf {k  >  Vl!jk  -  X)  , 


the  index  of  the  change  in  the  load  combination  process  which  is 
due  tot  the  nth’ change  in  the  constant  load  process.  Let 


v*> 


Ti  £  n 


1<it>  <Ti 


f 


the  number  of  that  occur  before  (discrete)  time  n  such 

-  x.  ,  =  k;  (lf.  ,(x)  =1  if  x  =  k  and  0  otherwise). 


that  x . 


For  each  integer  K 


%n  *  ui  1[0,K]<Ti  Vi  >  ' 


1  — 


the  maximum  load  combination  for  the  discrete  time  process  in 

the  time  interval  [0,n]  which  is  due  to  those  shock  loads 

and  constant  loads  which  occur  during  time  intervals  of  constancy 

for  the  discrete  time  constant  load  process  which  are  of  length 

(2) 

less  than  or  equal  to  K.  Finally,  let  ^  be  the  magnitude 
of  the  second  maximum  of  the  sequence 


(u.  1 


i  ~[0 


,Ki(Ti--  Vi1'  Tii.n> 


Hence 


Note  from  (3.2)  that 

i  Mn  • 


H  ■  (x)  i  P{M  <  x) 
K,n  K,n  - 


=  E[  n  [F*GK(X)J  n  ] 
k=0 


>  <  x) 


3.2.  Limiting  Results  for  Pareto-tail  Load  Distribution s 


F(x)  =  a.x  L (x)  ,  x  >  0  , 


G(x)  =  a2x  L  (x)  ,  x  >  0  , 


(3.5) 


where  L  is  a  slowly  varying  function. 

Next  we  will  find  the  limiting  distribution  of  as 

n  °°  ^or  fixed  K.  Then  we  will  show  that  the  limiting  distri¬ 
bution  of  is  close  to  that  of  for  K  large. 

PROPOSITION  (3.3).  Let  F  and  G  be  as  in  (3.5),  and  u 


be  such  that 


Then 


(u  )  aL(u)~i  as  n  . 
n  n  n 


(3.6) 


lim  H„  ( u  x)  =  exp{-x~a  p2  l  qk  [a,  +  ka~]>  =H.  (x) 
n-»-  °°  K,n  n  k=0  I  *  k 


Proof.  Note  that 


1  -  Gk(x)  =  kG(x)  G(x)] 


+  (k)  G(x) 2  G(x) k-2 


+  •••  +  [G(x5] 


~  ka2x  a  L(x) 


x  <» 
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■Thus  by  result  (8.13)  on  page  278  of  Feller  [1971] , 
1  -  F*G^  (x)  ~  (a1  -f  ka2)x”'a  L(jO  . 

..Hence  as  n  °° 

HK,n(unx)  =  EleXp{k^0  Ln(k)  £n  (unx)  }  1 


~  E[exp{  l  L  (k)  [1  —  F*GK(u_x)}] 
k=0  n  n 


~  ;E 


V(k) 

n 


x"a  [a1  +  ka2] 


Note  that 

L  (k)  2  k 

lim  — - —  =  p  q  ,  k  =  0,  1,...  . 

n 

The  result  now  follows  by  the  bounded  convergence  theorem. 


THEOREM  (3.4).  Under  the  assumptions  of  Proposition  (3.3) 

lim  P{Mn  1  unx^  =  exp{-x  “[pa^  +  qa2J)  . 
n  ■>  °° 

Proof.  Let 

K 

H*(x)  =  exp{  l  p2qkx”a [a.  +  ka?]} 

K  k=0  1  * 


and 


H (x)  »  lim  Uv(x) 
K  >  «  K 


=  exp{-x  ot[pa1  +  qa2]}‘. 


|p(Mn  <  v^x)  -  H(x)  | 


<  lP{Mn  <  .yd  -  HK(n(unx)l 


+  lHK_n(unx)  -  Hk(x)I  +  |HK(X)  -  H (x) i 


Note  that 


,P{Mn  -  unx}  "  HK,n(uhx)l  i  E U  -  exp{-An(K)}]  .  (3.7) 


where 


An(K)  = 


l  Ln<k)  An  F*Gk(u  x) 
k=K+l  n  n 


I  L  (k)  [1  -  F*Gk(u  x)  +  o(l  -  F*Gk(u  x))]. 
k=K+l  n  n  n 


Since 


I  i 


1  -  G  (x)  <  kG(x) 


1  -  F*Gk(x)  »  1-Gk(x)  +k  /  Gk_1(y)  F(x  -  y)  G(dy) 

0 

<  k  F*G(x) 


it  follows  that 


kn (K)  -  |  Ln'k>  U  s*G(*x)k  +  o(k  F*G(u  x)) ]  . 

iC— i\'T'X  ** 


Note  that,  for  0  <  e  <  3. 


F *  G(unx)  £F(unx(l-e))  +  G(unx(l-e))  +  F(unxe)  G(unxe)  . 
Therefore  for  fixed  6  >  0  there  is  an  Nn  such  that  for 


n  >  No 


p  “  L  (k) 

JL  VW  k  i  kJtl-V-  k  [*!*'“  +  a2x-“  +  «  1 


Further, 


[Vk)1  iEtJ1  ^k}'1!  *  *  »pV 


Hence  for  n  >  N_(6  ) 


EI*  l  L  <k)  k  F*5(u  x)] 


k?K+l 


lnx)  J  £  [a^^  +  a^Jx  +6]  J  p  q*k.  (3.8) 

k=K+l 


t-hoose  y  >  0.  Apply  Jensen's  inequality  to  (3.7): 


for  n  >  Nq  (5  ) 


E[1  -  expf-A^K)}]  £  l-exp{E[-A  (K)]}  0 


as  K  «.  Thus  for  sufficiently  large  n  and  K 


lP<Mn  <  V)  -  HK>n<V!> 


is  arbitrarily  small.  It  follows  that 


lim  P(M  £  xu  }  =  +  <3a2^ 

n 


(3,9) 


whereas,  if  the  load  process  were  one  of  independent  random  vari¬ 
ables,  denoted  by  with  distribution  pF(x)  +  qF*G(x),  the 

maximum  would  be  distributed  as  follows 


lim  P{M^  £  xun>  *  exp{-x*"ot[a1  +  qa2]}  • 

n  00 


(3.10) 


The  uR  is  the  same  in  each  case;  see  (3.6).  Comparison  of 
(3.9)  and  (3.10)  shows  that  the  expression  (3-10)  overestimates 
the  probability  of  exceeding  a  given  stress  level. 

O'Brien  [1974a]  obtained  the  above  result  in  the  case 
in  which  there  is  a  constant  load  (i.e.,  a2  =  0)  alone,  or 
only  a  shock  load  (i.e.,  a^  =  0). 


. . . . . . . . . 


3*3,  Results  for  Continuous  Time 

THEOREM  (3.5).  Let  F  and  G  be  as  in  (3.5)  .  Let  u(t) 

be  such  that  u(t)““  L( u(t))^  l/t  as  t  -►  *>.  Then  X+y. 

%  * 

lim  P{M(t)  <  xu(t)}  =  expf-x"’01  [Xa.  +  ya5] }  .  (3.9) 

t  -*- «  ~  x  4  . 

Proof.  Let  N(t)  be  the  number  of  changes  in  the  load  combina¬ 
tion  process,  (Z(u);  u  £  0),  in  the  time  interval  [Q,t).  The 
process  {N(t);t  >  0}  is  a  Poisson  process  with  rate 
Further,  • 

*<«  '*  «N(t)  •' 

The  result  now  follows  from  the  proofs  of  Proposition  (3.3) 
and  Theorem  (3.4)  ,,  the  strong  law  of  large  numbers  for  N(t), 
and  the  bounded  convergence  theorem. 

Slight  modifications  of  the  proofs  of  Proposition  (3.3) 
and  Theorems  (3.4)  and  (3.5)  imply  the  following  result. 

THEOREM. (3.6) ,  Let 

F  ( x)  =  x“aL(x)  and  G(x)  =  x_S  L(x) 

where  L  is  a  slowly  varying  function. 

a)  If  a  <  8  and  u(t)  is  such  that  u(t)  ”aL{u(t)  ]  ~  i/t  as  t  ■*> 

then  r  ,  , 

lim  P(M(t)  <_  xu(t) )  =  exp{-Xx  }.  (3.10) 

t  °° 
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b)  If  3  <  a  and  u(t)  is  such  that  u(t)  L[u(t)  1/t  as  ,  -t 
then 

lira  P(M(t)  £  xu(t)}  =  exp{-  yx"^  } . 
t  -*■  °° 

The  above  result  indicates  that  if  the  tail  of  the  dis¬ 
tribution  of  the  magnitude  of  a  constant  load  (respectively  shock 
load)  dominates  that  of  the  shock  (respectively  constant  load) , 
the  asymptotic  behavior  of  the  maximum  load  combination  is  the 
same  as  that  of  the  maximum  constant  load  (respectively  shock 
load)  by  itself.  On  the  other  hand,  Theorem  (3.5)  indicates 
that  if  tails  of  the  distributions  of  the  magnitudes  of  the 
constant  and  shock  loads  are  comparable,  then  the  asymptotic 
behavior  of  the  maxima  of  the  load  combination  process  depends 
on  both  the  constant  and  shock  load  processes . 

3.4 .  Limiting  Results  for  the  Joint  Distribution  of  the  First 
and  Second  Maxima 
(2) 

Let  Mn  denote  the  magnitude  of  the  second  maxima  of 
(2^;  k  £n}.  Assume  F  and  G  are  as  in  (3.5).  Let 

H (x)  =  exp{-x_a[pa1  +  qa2J} 


as  before.  Techniques  similar  to  those  used  in  the  proofs  of 
the  results  in  Subsection  (3.2)  yield  the  following  result 


THEOREM'  (3.7)  .  Let'  u‘  be  as1  in  (3.6)  .  Then 

n 

lim  P(M^2)  <  yun,  <  xun}  =  B(y,x) 


where*  for  x  £  y, 


B(y,x)  =  H  (x) 


and'  for  y  <  x 


B(y,x)  =  H(y)  {1  +  [p2a1  +  qa2]  [y~a  -  x"a]}  . 


A  result  of  Welsch  [1972]  implies  that  the  limiting 
distribution  B  must  be  of  the  following  form 


B(y,x)  = 


•  H  (x) 


if  y  >  x. 


H(y){l-g[(£n  H(x)/£n  H(y))]£n  H(y)}  ,  if  y  <  x 


where  g(s-)  ,0<s<^l  is  a  concave,  nonincreasing  function  which 
-satisfies  g(0)(l  -  s)  <_  g(s)  <_  1  -  s.  In  our  case 


g(s) 


"  P2ax  +  qa2“ 
pa,  +  qa9 


(1  -  s)  . 


(3.12) 


Let  M^  (t)  denote  the  magnitude  of  the  second  maximum 
to  occur  in  [0, t]  for  the  continuous  time  load  combination  . 
process  (Z(s);  s  J>  0}  .  Arguments  similar  to  those  in  Section 
3.3  can  be  used  to  obtain  the  following  result. 

THEOREM  (3.8).  Let  F  and  G  be  as  in  (3.5)  and  u(t)  be 
such  that  uft)”01  L(u(t))  ~  1/t  as  t  +  «.  For  x  >  y 


lim  P(M(t)  £xu(t),  KU)  (t)  £  yu(t)  } 
t  °° 


=  exp{-y  atAa^  +  uajl)  1  + 


for  x  £  y 


a.  +  ya,  (y  -  x 


A+^I  “1  '  ^  2 


■“)]  ; 


lim  P{M(t)  £  xu(t) ,  M(2) (t)  £  yu(t) ) 


t  ■>  °° 


=  exp{-x  a[Aa^  +  ya2J}» 
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